Topological quantum computing with a noisy network and percent-level error rates 



Naomi H. Nickerson 1 , Ying Li 2 , Simon C. Benjamin 2 ' 3 

1 Department of Physics, Imperial College London, Prince Consort Road, London SW7 2AZ, UK 
2 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543 and 
3 Department of Materials, University of Oxford, Parks Road, Oxford 0X1 3PH, UK. 
Please address correspondence to: s.benjamin@qubit.org 

A quantum computer could be built by networking together many simple processor cells. The 
difficulty is that realistic quantum links are very error prone. A solution is for cells to repeatedly 
communicate with each other and so 'purify' any imperfections; however prior studies suggest that 
the cells themselves must then have extremely low internal error rates. Here we describe a method 
by which even error-prone cells can perform purification: groups of cells generate shared resource 
states, which then enable stabilization of topologically encoded data. Given a realistically noisy 
network (> 10% error rate) we find that intra-cell error rates for initialisation, state manipulation 
and measurement can simultaneously exceed 0.82% before the protocol fails. 



Topological codes are an elegant and practical method 
for representing information in a quantum computer. 
The units of information, or logical qubits, are encoded 
as collective states in an array of physical qubits pQ . We 
employ stabilizers, i.e. measurements on nearby groups 
of physical qubits, to detect errors as they arise. By a 
suitable of choice stabilizers we can even generate logical 
operations between the encoded qubits. The act of ap- 
plying stabilizers over the array thus constitutes a kind of 
'pulse' for the computer - it is a fundamental repeating 
cycle and all higher functions are built upon it. 

If the stabilizer operations themselves are perfect, then 
a topological code can be remarkably robust to noise; 
for example ~ 19% of physical qubits can suffer errors 
before the encoded information is corrupt [2]. However 
here we are concerned with the case that the stabilization 
process itself can suffer errors, through imperfections in 
steps such as initialising ancillas, performing gate oper- 
ations and measurements. Stabilizers will still "do more 
good than harm" , and extend the lifetime of encoded in- 
formation, provided that such errors are sufficiently infre- 
quent. Estimates of the threshold error rate are around 
0.75% - 1.4% depending on the model [M5]. 

These numbers are applicable to a monolithic architec- 
ture (Fig. [I] upper), i.e. a large grid of qubits with each 
connected to its neighbours such that two-qubit opera- 
tions can be performed with high fidelity 5j. Building 
such a structure is a plausible goal for some systems, no- 
tably specific superconducting devices [Sj. However it is 
less appropriate for other nascent quantum technologies. 
For example an individual nitrogen-vacancy (NV) cen- 
tre in diamond, with its electron spin and small number 
of associated nuclear spins, naturally forms an isolated 
cluster [7j. Similarly a small ion trap holding a modest 
number of ions is a self-contained unit [5]. We use the 
term 'cell' for these convenient small groupings of qubits. 
Cells could be linked up to form a networked array (Fig. [I] 
lower), but we should expect links to have relatively poor 
fidelity. Links may be photonic, in which case a key issue 
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FIG. 1: Quantum architectures. Upper figure follows the 
scheme of Ref. [5]. These schematics are examples of very 
broad classes of systems; e.g. the network may actually be 
composed of NV centre 'cells' linked by spin chain 'wires'. 



is photon loss, and instabilities in path lengths or inter- 
action strengths can also degrade the channel. An alter- 
native approach in the solid state is to employ a chain 
of spins to 'wire' together different cells [S]. Here again 
it may be infeasible to achieve error rates below the 1% 
level, given realistic decoherence affecting the spins [lOj . 

Fortunately the problem of poor inter-cell links is not 
insurmountable: cells can interact repeatedly and purify 
the results to remove errors [TTJ [T3] . Thus one can realise 
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a 'noisy network' (NN) paradigm, which is also called dis- 
tributed |13l I14j or modular |15) . However prior analyses 
of this approach have indicated a serious drawback. The 
need to perform purification means a cell's own internal 
error rate must be very low; of the 0.1% order [TTMl5| . 
This contrasts poorly with the 1% general error thresh- 
old for the monolithic case. Given that a real machine 
should operate well below threshold, such a demand on 
the intra-cell precision may be prohibitive. 

Here we describe a new NN protocol which achieves a 
threshold for the intra-cell operations which is far higher, 
and comparable to current estimates for tolerable noise in 
monolithic architectures. We introduce an optimised net- 
work protocol enabling our cells to collectively represent 
logical qubits according to the 2D toric code (one of the 
simplest and most robust topological surface codes Q], 
directly related to the approach used in cluster state ar- 
chitectures [3l fT3Hl5l lT7]). This code is stabilized by re- 
peatedly measuring the parity of groups of four qubits in 
either the Z or the X basis. We implement such measure- 
ments by directly generating a GHZ state shared across 
ancilla qubits in four cells, before consuming the resource 
in a single step to jointly stabilize one "data qubit" in 
each cell. This procedure is in the spirit of the bandaid 
protocol |16j and contrasts with the standard approach 
of performing a sequence of two-qubit gates between the 
four data qubits and a fifth auxiliary qubit (which, in the 
network picture, would require its own cell). By directly 
generating the GHZ state over the network we remove 
the need for the auxiliary unit and more importantly we 
can considerably reduce the accumulation of errors. 

The protocol specified in Fig. [2] requires a total of four 
qubits within each cell of the network. One qubit in each 
cell is designated the data qubit. The other three are 
ancillas that will be initialised, processed and measured 
during each stabilizer evaluation process. This proves 
sufficient to handle network error rates at the 10% level 
while still achieving a near-percent level threshold for 
intra-cell error rates. We note that if a further ancilla 
qubit is available in each cell then one can tolerate still 
higher levels of network noise, or alternatively one can 
boost the intra-cell threshold even closer to 1% (please 
see appended Sup. Mat. document). 

We use the network channel between two cells to create 
shared, noisy Bell pairs in the Werner form 
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where |$+) = (|00) + \11))/V2. Throughout this paper 
we choose the network error probability to be p n = 0.1. 
This 'raw' Bell generation is the sole operation that oc- 
curs over the network. We additionally require only three 
intra-cell operations: controlled- Z (i.e. the two-qubit 
phase gate), controlled- AT (i.e. the control-NOT gate) 
and single qubit measurement in the A-basis. The two 
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FIG. 2: The three phases for generating and consuming a 
4-cell GHZ state among the ancilla qubits (blue-to-green) in 
order to perform a stabilizer measurement on the data qubits 
(purple). Phase (1): use purification to create a high quality 
Bell pair shared between cell A and cell B, while in parallel 
doing the same thing with cells C and D. (2) Fusion opera- 
tions, A to C and B to D, create a high fidelity GHZ state. 
(3) Finally use the GHZ to perform a one-step stabilizer op- 
eration. The parity of the four measured classical bits is also 
the parity of the stabilizer operation we have performed on 
the data qubits. Two dashed regions indicate operations that 
are part of the STRINGENT protocol; omitting them yields 
the EXPEDIENT alternative. 



qubit gates are modelled in the standard way: an ideal 
operation followed by the trace-preserving noise process 

N( P ) = (i - Pg )p +^j2( A ® B )p( A ® B ) f 

10 A,B 

where operator A € {1, a x , o- y ,a z } acts on the first qubit, 
and similarly B acts on the second qubit, but we exclude 
the case 1 <g> 1 from the sum. Meanwhile noisy measure- 
ment is modelled by perfect measurement preceded by 
inversion of the state with the probability p m . Phases 1 
and 2 in Fig. [2] are posts elective: whenever we measure a 
qubit, one outcome indicates 'continue' and the other in- 
dicates that we must recreate the corresponding resource 
(the desired outcome is that which we would obtain, were 
the noise parameters set to zero). We note that in the 
circuits developed here we adopt and extend the "double 
selection" concept introduced by Fujii and Yamamoto in 
Ref. [18] . The Figure shows two variants of our protocol. 
EXPEDIENT minimises the number of raw pairs required 
and thus the overall time requirement. On average it re- 
quires fewer than twice the number of pairs consumed in 
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the ideal case that all measurements "succeed". Alter- 
natively the STRINGENT variant imposes higher fidelity 
standards at the cost of a longer procedure. The perfor- 
mance of both protocols is found presently. 

Having specified the stabilizer protocols, we must de- 
termine their real effect given the various error rates p n , 
p g and p m . It is convenient derive a single superopera- 
tor describing the action of the entire stabilizer protocol 
on the four data qubits. In the following we write M to 
stand for the "odd" or "even" reported outcome of the 
stabilizer protocol (i.e. the parity of the four measured 
qubits in Phase 3 of Fig. [2|. If density matrix p represents 
the state of all data qubits prior to the evaluation of our 
stabilizer, then measurement outcome M and the cor- 
responding state S M {p) = P M (p)/Tr[P M (p)} will occur 
with probability Tr[P M (p)], for some projective opera- 
tor P M {-). Now stabilising the toric code involves two 
types of operation, the X and the Z— stabilizers. Sup- 
pose we are performing a Z stabilizer. If our protocol 
could act perfectly, then for example the even projec- 
tor would be -Pide™(p) = 12 l*)(*|p|j)0'|) where the sum 
is over all states \i), \j) with definite even parity in the 
Z-basis: |0000), |0011), etc. Analogously the ideal odd 
Z-stabilizer sums over states of definite odd parity, and 
meanwhile for X stabilizers the ideal projectors refer to 
states of definite parity in the A-basis. In reality our 
imperfect operations result in projectors of the form 

P?L(P) = E aeE e P^{p)El + bfE e P™ eal (p)El 

e 

with E e = (ABCD) e and {A, B, C, D} £ {1, <r x , a y ,a z }. 

Here the four operators making up E are understood to 
act on data qubits 1 to 4 respectively, and index e runs 
over all their combinations. The symbol M represents 
the compliment of M, i.e. "odd" for M ="even" and 
vice versa. We see that this real projector is made up 
of a mix of the correct and incorrect projectors together 
with possible Pauli errors; the various weights a and b 
capture their relative significance. 

Given the underlying error rates p n , p g and p m we 
can employ the Choi-Jamiolkowski isomorphism to find 
the corresponding weights {a, 6} in our superoperator. 
In the Supplementary Material we give examples. We 
find that the largest contributor to P^aiiP) IS * ne re ~ 
ported parity projection, i.e. for M ="even", the largest 
of all weights {a, b} is that associated with E = (1111) 
and -PidoaK/ 5 )- The next largest term will be the pure 
"wrong" projection, i.e. the combination of E = (1111) 
and ^dcai(/°)- This form of error is relatively easy for the 
toric code to handle, and we have deliberately favoured 
it over other error types by minimising the o x and o~ y 
errors, rather than cr z , in our stabilizer protocol. The re- 
maining terms correspond to Pauli errors in combination 
with either the correct or the incorrect projectors; for 
example E — (o x 111) is an erroneous flip on data qubit 
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FIG. 3: Scheduling stabilizer operations. The right side 
graphic shows the standard arrangement of one complete sta- 
bilizer cycle, involving Z and X projectors (square symbols 
indicate that the four surrounding data qubits are to be sta- 
bilised). Because a given cell can only be involved in gen- 
erating one GHZ resource at a given time, each of these two 
stabilizer types must be broken into two subsets; see main fig- 
ure. Fortunately in our stabilizer superoperator S(p) we can 
commute projectors and errors so as to expel errors from the 
intervening time between subsets, so allowing them to merge. 



I. Such single-qubit errors are more probable than two- 
qubit errors, and the occurrence of three or four errors is 
highly improbable. 

Having determined the stabilizer superoperators 
^roai(') f° r both Z-type and X-type stabilizers, together 
with a suitable scheduling scheme as shown in Fig. |3j 
we proceed to determine thresholds by intensive numer- 
ical modelling. Our model introduces errors randomly 
but with precisely the correlation rates indicated by the 
weights {a, b}. We record the (noisy) stabilizer outcomes 
and subsequently employ Edmonds' minimum weight 
perfect matching algorithm [TH1 HO] to pair and resolve 
stabilizer flips in the standard way (see e.g. Ref. [S]). 
This technique allows us to establish the threshold for 
successful protection of quantum information by simu- 
lating networks of various sizes. If an increase of the 
network size allows us to protect a unit of quantum infor- 
mation more successfully, then we are below threshold. 
Conversely if increasing the network size makes things 
worse, then in effect the stabilizers are introducing more 
noise than they remove and we are above threshold. The 
results of a large number of such simulations are shown 
in Fig. [5] We see thresholds for EXPEDIENT and for 
STRINGENT of over 0.6% and 0.77% respectively. The 
third graph shows the performance of a five-qubit-per-cell 
variation of EXPEDIENT which we describe in more de- 
tail in the Supplementary Material. It achieves a thresh- 
old in excess of 0.82% by employing an additional filter 
such that most stabilizers are improved while a known 
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FIG. 4: Performance of the EXPEDIENT and STRINGENT 
network protocols, given a toric code with n rows x n columns 
(2n 2 qubits) of data qubits. A given numerical experiment is 
a simulation of 100 complete stabilizer cycles on an initially 
perfect array, after which we attempt to decode a Z measure- 
ment of the stored qubit. The result is either a success or a 
failure; for each data point we perform at least 10,000 exper- 
iments to determine the fail probability, and reciprocate this 
to infer an expected time to failure. Network error rates are 
10% in all cases; we set intra-cell gate and measurement error 
rates equal, p m — p g , and plot this on the horizontal axes. 
For low error rates the system's performance improves with 
increasing array size. As the error rates pass the threshold this 
property fails. Insets: typical final states of the toroidal ar- 
ray after error correction. Yellow squares are flipped qubits, 
green squares indicate the pattern of Z-stabilizcrs. Closed 
loops are successful error corrections; while both arrays are 
therefore successfully corrected, it is visually apparent that 
the above-threshold case is liable to long paths. 



minority become more noisy, and this classical informa- 
tion is fed into an enhanced Edmonds' algorithm. 

The Edmonds' algorithm approach employed here has 
been chosen because it has been well studied in the con- 
text of noisy stabilizers; in the Supplementary Material 
we apply our same model to monolithic architecture, ob- 
taining a threshold in the region of 0.9%, generally con- 
sistent with prior studies [5]. We note that for the ide- 
alised case of perfect (noiseless) stabilizers other algo- 
rithms have been developed which may offer advantages 



over Edmonds' [31[3TJ[22]. Such approaches might even- 
tually translate to higher thresholds for noisy stabilizers 
(both for monolithic and noisy network cases). 

We have yet to consider general decoherence caused 
by the external environment during a stabilizer protocol. 
Fortunately many of the systems most relevant to the 
NN approach have excellent low-noise "memory qubits" 
available. For example in NV centres at room tempera- 
ture, nuclear spins can retain coherence for the order of a 
second; and for impurities in silicon the record is several 
minutes 23, 24]. We would naturally use such spins for 
our data qubits and for the innermost ancilla qubits, i.e. 
those bearing the GHZ-state in Fig. [2] Then given that 
a conditional gate on a nuclear spin may take of order 
10/ns, we estimate the probability of an environmentally 
induced error during a full stabilizer protocol to be at 
least an order of magnitude below the error rates due to 
p g and p m (see appended Sup. Mat. document). 

In conclusion, we have described an approach to 'noisy 
network' quantum computing, where many simple mod- 
ules or cells are connected with error prone links. We 
show that relatively high rates of error within each cell 
can be tolerated, thus closing the gap between error tol- 
erance in networks versus monolithic architectures. We 
hope this result will encourage the several emerging tech- 
nologies for which networks are the natural route to scal- 
ability; these include ion traps and NV centres (linked 
either optically or via spin chains). 

We thank Ben Brown and Earl Campbell for numer- 
ous helpful conversations. This work was supported by 
the National Research Foundation and Ministry of Edu- 
cation, Singapore. 
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Supplementary Material 



COMPARISON TO THE MONOLITHIC CASE 

In order to make a clear comparison with the threshold 
that can be achieved with the monolithic architecture, we 
applied our same superoperator description and numeri- 
cal simulation to this case. (Of course there is no meaning 
for p n , the network error rate, in a monolithic architec- 
ture). The figure shows the circuit we used for the case 
of a Z stabilizer. Note that this approach requires ini- 
tialisation of the single shared auxiliary qubit; we used 
the value of p rn as the fidelity of this initialisation. 

One sees that the threshold is between 0.9% and 0.95% 
and is therefore appreciably higher than that obtained by 
STRINGENT whilst in the same "ball park" . As noted 
below, the addition of another ancilla per cell can further 
close this gap. 
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FIG. 5: Performance of the monolithic architecture, to be 
compared with the graphs in Fig. 4 of the main paper. 



TWIRLING 



The protocols given in the main paper lead to slight 
irregularities between the weights associated with given 
errors occurring on different specific data qubits; conse- 
quently in our superpoperator the weight associated with 
ZZU might be slightly higher or lower than the weight 
associated with tZZt . This is because the protocol per- 
forms A — B, C — D pairings in Phase 1, and then A — C, 
B — D pairing in Phase 2, thus the errors that survive 
this purification process will not be equally distributed 
under rotation of cell labels. While these irregularities 
have no significance for our threshold calculations, they 
do cause the weightings to have a spurious complexity. 
Therefore in our analysis we append onto the procotols of 
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Fig. 2 a additional twirling operation which randomly ap- 
plies swap operations between cells so as to 'smooth' the 
weightings. (Physically this is equivalent to the rather 
perverse act of programming one's system with a range 
of possible protocols, identical except for permutation of 
the cell labels A to D, and then applying one at random 
without retaining a record of the choice.) 

EXAMPLE SUPEROPERATORS 

In the main paper we wrote the following to represent 
the overall effect of the stabilizer protocol on four data 
qubits initially in state p, 

PriUp) = E a * E ePifeAP)4 + bf E e P^p)E^ e with 
e 

E e = (ABCD) e and {A, B, C, D} e {1, a x , a y ,a z }. 

Here M stands for the reported outcome, "odd" or 
"even", and the -P^aiO) represents the perfect parity 
projector in either the Z or X basis depending on which 
class of stabilizer one is performing. The four operators 
making up E are understood to act on data qubits 1 to 
4 respectively, and index e runs over all their combina- 
tions. The symbol M represents the compliment of M, 
i.e. "odd" for M = "even" and vice versa. 

We noted that this real projector is made up of a mix 
of the 'correct' and 'incorrect' ideal projectors together 
with possible Pauli errors; the various weights a and b 
capture their relative significance. Here we give some 
examples of those weightings. The full list has 256 terms 
but they rapidly fall in magnitude so that the latter half 
are extremely small; here we list the terms corresponding 
to all single- and two-qubit errors. For the examples we 
will give, in fact the superscript M can be omitted from 
the a and b weights, i.e. the same set of weights a, b 
apply in P^™ as in P°^. In other words the weight 
of the 'correct' projector P^ ai is the same in each, as are 
the weights of all the various error combinations. We will 
present the numbers as follows 

{{At,A z , Ax,A y ,A Z z,A xx ,Ayy,Axz, A yz , A xy }, 
{.Bi, E>z, Bx,E>Y-, Bzz,Bxx,Byy,Bxz, Byz,Bxy}} 

Here the weights A are associated with the 'correct' pro- 
jector -P^cai an d the B weights are associated with the 
'incorrect' projector P^^- The subscripts refer to the 
Pauli errors in the corresponding term, so that the sub- 
script 1 indicates no errors, while subscript X means 
"a ax on one of the four data qubits" , and YZ means 
"a oy on one of the qubits and a oz on another". 
These capital letters A and B are therefore sum of one 
or more of the weights a, b in Eqn. [T] for example 
Ax = a llix + aiixi + a 1X tt + axiti- As noted in 
the previous section, the twirling technique means that 



the individual a or b terms in such a sum in fact all have 
the same value. 

For the case of EXPEDIENT operating with p n = 0.1 
and p m = p g — 0.006 we find projectors in the Z basis 
have following weights: 

{{0.9117, 0.00681, 0.00314, 0.00314, 0.000182, 0.00000758, 
0.00000758, 0.0000336, 0.0000336, 0.0000152}, 
{0.0617, 0.00674, 0.00314, 0.00314, 0.000127, 0.00000758, 
0.00000758, 0.0000336, 0.0000336, 0.0000152}}. 

(Meanwhile projectors in the X basis have the same num- 
bers if we exchange subscripts If>Zin the listing or- 
der used above.) Note that many of the corresponding 
A and B terms are identical due to symmetries in the 
algorithm (including the twirling). We see that there is a 
91% chance of a pure 'correct' projection, a 6.2% chance 
of a pure 'wrong' projection, a 0.68% chance of a correct 
projection followed by a single az error, and so on. For 
the case of STRINGENT operating with p n = 0.1 and 
Pm = Pg = 0.0075 we find, 

{{0.928, 0.00675, 0.00391, 0.00391, 0.000187, 0.00001182, 
0.00001182, 0.0000414, 0.0000414, 0.0000236}, 

{0.0424, 0.00665, 0.00391, 0.00391, 0.0000849, 0.0000118, 
0.00001182, 0.0000414, 0.0000414, 0.0000236}}. 

Finally for the comparison case of the monolithic circuit 
where there is no network noise to consider, setting p. m — 
p g = 0.09 we find 

{{0.951, 0.00470, 0.00352, 0.00352, 0.00119, 0.0000128, 
0.0000128, 0.00121, 0.00121, 0.0000256}, 
{0.0178, 0.00470, 0.00352, 0.00352, 0.00119, 0.0000128, 
0.0000128, 0.00121, 0.00121, 0.0000256}}. 

FIVE QUBITS PER CELL 

The protocols in the main paper use a total of four 
qubits per cell of the network (one data qubit, three an- 
cillas). If one adds further ancilla(s) then the perfor- 
mance improves. While it is beyond the scope of the 
present paper to seek optimal protocols that exploit five 
(or more) qubits, we can easily modify our approaches 
to make some use of the additional resource. For exam- 
ple, wherever our original 4-qubit STRINGENT protocol 
calls for a "raw" Bell pair to be created on a given an- 
cilla pair we can instead insert a small circuit that creates 
raw pairs using the additional ancillas and purifies them. 
Similarly one can replace "single selection" purification 
(two tiers of ancilla yield an improved pair on the higher 
tier) with "double selection" (three tiers of ancillas yield 
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Round two Bell pair production 


6 


0.7506 


1 


3 


Round one single selection 


4 


0.8619 


3 


4 


Round two single selection 


3 


0.8550 


3 


5 


Make GHZ 


2 


0.8651 


1 


6 


Round one single selection 


4 


0.8619 


6 


7 


Round two single selection 


3 


0.8550 


7 


8 


Check GHZ 


2 


0.8654 


1 


9 


Measure stabilizer 


2 







TABLE I: Listing of probabilities used in modelling the number of steps required for EXPEDIENT. 



a significantly improved pair on the highest tier) . Adopt- 
ing these rather naive modifications we immediately find 
that the network noise is tolerable at the p n — 0.2 level 
(rather than p n = 0.1) for the same threshold of 0.77% 
local errors. 

Alternatively we can hold the network errors constant 
and raise the intra-cell error rate. Data for this case are 
shown in the third panel of Fig. 4 in the main paper, 
labelled "STRINGENT+" . We have introduced one fur- 
ther enhancement: We replace the usual Phase 3 of the 
STRINGENT protocol i.e. the steps where the GHZ re- 
source would simply be coupled to the data qubits and 
then measured out. Instead we filter the GHZ after cou- 
pling it to the data qubits, and if this filter fails we abort 
the protocol by measuring the GHZ in the Z-basis (this 
therefore requires the addition of Z-basis measurement 
to our set of allowed primitive operations). If we have 
aborted then we then perform a whole new round of sta- 
bilization, this time without the filter so that the protocol 
will certainly complete. This procedure leaves us know- 
ing the stabilizer outcome and some additional classical 
information about exactly what steps occurred. Specif- 
ically there are 3 distinct cases: (a) The filter was suc- 
cessfully passed; this is the best case and results in lower 
error rates on the data qubits it is about 92% of cases 
in for the parameter range we consider, (b) The filter 
failed, but on measuring the GHZ in the Z-basis it was 
found to be in a correct state (e.g. 0000 or 1111). This 
occurs about 4% of the time. It is another "good" case 
in that there is a low chance of errors having reached 
the data qubits, so that the second round can perform 
normally, (c) The filter failed, and on measuring on the 
GHZ in the Z-basis it was found to be in an incorrect 
state (e.g. 0001). This occurs about 4% of the time and 
it is the "bad" case; in this event there is very likely to 
be an error on the data qubits. Now if we make no use of 
the classical information and merely "forget it" then the 
net effect of this protocol is to make things worse versus 
a simple one-round use of STRINGENT - this is not sur- 
prising since the overall risk of an error is not reduced (it 
is increased slightly due to the extra steps). However if 



we modify our Edmonds matching algorithm to use the 
classical information, specifically to favour paths that are 
consistent with errors occurring where "bad" stabilizers 
took place, then the threshold improves (this is the case 
shown in the third panel of Fig. 4). In effect we trade 
a small amount of increased error risk for a significant 
amount of classical knowledge: the 4% of "bad" stabi- 
lizers account for about half of all errors entering the 
system. 



MEMORY ERRORS 

Since our protocols are considerably longer than the 
circuit required for the monolithic architecture, and in- 
deed our approach is post-selective and therefore of un- 
certain duration, it is important to assess the potential 
impact of memory errors. In the main paper we assert 
that memory errors should have negligible impact, if our 
cells can employ qubits at least as good as those demon- 
strated in the Science papers of Maurer et al and Steger 
et al. To substantiate this we need to estimate the dura- 
tion of the protocols; here we do so for EXPEDIENT since 
it is the protocol designed for use when memory errors 
are an issue. In the following we assign one "time step" 
to any elementary operation in our protocol, whether 
a gate, a raw Bell creation or a measurement. We as- 
sume that the operations within a cell must be strictly 
sequential. Because we are evaluating an entire 'sheet' of 
stabilizers in parallel across the array (see Fig. 3 of the 
main paper), we need to be concerned not merely with 
the average time that a stabilizer protocol might take, 
but rather with the the time required for a given target 
proportion of all stabilizers to succeed. Our simulations 
indicate that if we wait for 99% of stabilizers to evaluate, 
and abandon the remaining 1%, then there is negligible 
impact on the threshold. (A stabilizer that fails to report 
is less damaging than a stabilizer that performs a 'wrong' 
projection; moreover such stabilizers will not introduce 
errors since the ancilla GHZ state is never created and 
coupled to the data qubits.) Therefore we simply take 
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the expected time for 99% of stabilizer protocols to com- 
plete as the characteristic time for a 'sheet' of stabilizers 
to be found. Obviously there is scope for far more sophis- 
ticated approaches which minimise waiting time, but this 
naive strategy suffices to give us a bound. 

The process of building a GHZ state can be separated 
into distinct sections, each of which is terminated by a 
measurement, which, if it results in the 'wrong' outcome, 
will reset the process to an earlier stage. Each of these 
measurement outcomes has an associated probability of 
success, pl, where the index L denotes the level. For 
the EXPEDIENT protocol operating at error rates of 
P g = Pm — 0.6% and p n = 10%, these levels and their 
probabilities are given in Table [TJ It should be noted that 
levels 1 and 2 utilise two cells (and 3,4 and 6,7) are run 
twice in parallel accross four cells, and consequently the 
longer of the two times will determine when the process 
may proceed to the next level. 

The GHZ production process is probabalistic, with a 
minimum duration of — 33, which occurs with 

a probability JliPi = 0.2242. Using the parameters 
in the table the process was simulated; 100,000 samples 



were generated and used to estimate the parameters of 
the resulting distribution. This found the expectation 
time per stabilizer to be 68.2 time steps, while 50% of 
operations were completed after 57 time steps, 95% after 
138 time steps, 99% after 195 time steps and 99.9% after 
278 time steps. 

Taking the value 195 steps, we find that 24 of these 
steps are manipulations of the memory qubits, i.e. a 
data qubit or a GHZ-bearing tier of the ancillas, which 
we are now taking to be nuclear spins. We can assume 
that conditional gates on these qubits are by far the slow- 
est, at perhaps 10/zs (see the Steger reference in the main 
paper). We therefore take 500/is as a characteristic time 
for a 99% subset of a complete 'sheet' of stabilizers to 
evaluate. Maurer et al reported that their best NV mem- 
ory qubit had lifetimes of about 2 seconds (note the sili- 
con impurity qubits were reported to live far longer, ~ 3 
mins) . If a proportion 1/e survive 2 seconds, we infer an 
error rate about 40% per second, or about 0.025% over 
500/is. This can indeed be neglected in comparison to 
the error rates from the active processes, i.e. the gates 
and measurements. 



